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Abstract. A vital ingredient of DNA molecule named adenine may be produced 
by successive addition of HCN during molecular cloud collapse and star formation. 
We compute its abundance in a collapsing cloud as a function of the reaction rate 
and show that in much of the circumstances the resulting amount may be sufficient 
to contaminate planets, comets and meteorites. We introduce a /-parameter which 
may be used to study the abundance where radiative association takes place. 



Keywords: : Molecular Clouds, Star Formation, Organic Compounds, Biomolecules 
PACS Nos. : 98.38.Dq, 97.10.Bt, 61.66.Hq, 87.15.R 

: Accepted for Publication as RAPID COMMUNICATIONS for April 1, 
2000, issue of Indian Journal of Physics. 1. Introduction 

In recent papers, Chakrabarti [1], and Chakrabarti and Chakrabarti ([2], 
hereafter referred to as Paper I) explored the possibility of the formation 
of biomolecules in star formation region using gas-phase chemistry. Their 
conclusion was that even in frigid condition in interstellar matter some of 
the simplest amino acids such as glycine, alanine etc. could be produced 
even before the formation of stars and planets. Paper I also showed that 
with a choice of reaction rate constant 10 -10 , significant adenine may also 
be produced. Some preliminary results of amino acids are in [1]. 
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Adenine is a simply produced vital component of the DNA molecule and 
its significant production may point to an important clue into the problem 
of origin of life on planets like ours. Because of this, it is essential to carry 
out careful analysis on the reaction rates during the adenine formation. 
In Paper I, we used the formation of adenine by successive addition of 
HCN [3] by using an 'average' rate. In the normal circumstances, in gas- 
phase reaction HCN + HCN — ► H2C2N2 rate would be small, since they 
must combine by radiative association, i.e., they must radiate a photon 
when combined together. This is a slow process and the probability of 
photon emission could be 1 in a few thousand to a few million (T. Millar, 
private communication). However as the size of the molecule gets bigger, the 
process becomes faster. Thus, it is likely that for a large enough molecule, 
the radiative association may take place at every collision and at this stage, 
the collisional rate may be used. One possibility is therefore to assume that 
after every addition of HCN, the reaction rate goes up by a factor of / (/ 
may be anywhere from 1 to 100 or more). Hence one may imagine that at the 
early stages, HCN '+ HCN — ► H2C2N2 forms with a reaction rate of 10 -16 , 
but for HCN + H 2 C 2 N 2 the rate becomes / x 10" 16 , for HCN + H 3 C 3 N 3 
the rate becomes f 2 x 10~ 16 and so on. It would be therefore interest to 
learn, whether significant adenine is formed and it is detectable when the 
radiative association process is taken into account. In the present paper, 
we do just that. It is possible that more favorable reactions take place on 
ice, but in view of little known reaction rates of ice chemistry we believe 
that the best we could do is to study the formation of these important 
molecules as a function of two parameters, namely, a Ad and /. It is quite 
possible that a suitable / parameter we suggested above would take care 
of the ice-chemistry reaction rates as well. It is still possible that such an 
/ may actually be determined by actual detection of molecules in space. 
Similarly, constancy of / is an assumption of our model. In reality it could 
vary with the size of the molecules. 



So far, there has been controversy whether glycine has been observed 
in interstellar matter. Miao et al. [4] tentatively detected glycine in the 
massive star forming region Sgr B2(N) though this was later challenged by 
Combes et al. [5] who suggested that with the sensitivity of the detector 
taken into account, the lines were really at the confusion limit and positive 
identification would require more sensitive instruments. It is not known 
if any attempts were made to detect adenine lines, however there may 
have been detection of adenine in meteoritic samples (M. Bernstein, private 
communication) . 
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2. Reaction Network and Hydrodynamic Model 

We choose the same reaction network as in Paper I which we again present 
here for the sake of completeness. We take the UMIST database (Millar, 
Farquhar & Willacy [6]) as our basis of chemical reactants and reactions, 
but added several new reactions such as synthesis of amino acids (alanine 
and glycine), hydroxy-acids (glycolic and lactic acids), DNA base (adenine, 
[3]), urea synthesis etc. These new reactions make the total number of 
species to be 422. The rate co-efficients of these additional reactions are 
difficult to find, especially in the environs of a molecular cloud. To use 
UMIST database, the rate constant for a two body reaction is written as 
[6], 

k = a(T/300) /3 exp(- 7 /T) cm 3 s" 1 (1) 

where, a, (3 and 7 are constants and T is the temperature. Amino acid syn- 
thesis rate was estimated from Fig. 8 of Schulte & Shock [7]. Urea synthesis 
rate is kept comparable to the rates given in UMIST table. The rate con- 
stants were taken to be a = 10~ 10 , j3 = 7 = for each two-body reactions. 
In Paper I, the rate constants for adenine synthesis was chosen to be similar 
to other two body reactions [ctAd = 1CP 10 (3 = 7 = for each HCN addition 
in the chain HCN -» (on) -» CH(NH)CN -> (a 2 ) -> NH 2 CH(CN) 2 -> 
(a 3 ) -» NH 2 (CN)C = C{CN)NH 2 -> (a 4 ) -» H 5 C 5 N 5 (adenine)]. In 
the present paper we run the same simulation with on |j=i, 2,3,4 = a Ad = 
10~ 12 , 10~ 14 , 10 -16 as well and in addition, consider the possibility that 
on = f 1-1 a Ad as discussed in the introduction. In this notation, Paper I, 
represents the case with / = 1. 

Initial composition of the cloud before the simulation begins is kept to 
be the same as in [6], and formation of H 2 is included using the grain-surface 
reaction with rates as in [6]. The initial mass fractions are taken to be the 
same as in [6] (but appropriate convertion), i.e., H:He:C:N:0:Na:Mg:Si:P:S:Cl:Fe 
= 0.64:0. 35897:5. 6x 10" 4 :1.9x 10" 4 :1.81x 10" 3 :2.96x 10" 8 : 4.63x 10~ 8 :5.4x 
10" 8 :5.79 x 10~ 8 :4.12 x 10~ 7 :9 x 10~ 8 :1.08 x 10~ 8 . 

The hydrodynamic model is kept same as that in Paper I. We choose 
the initial size of the molecular cloud to be tq = 3 x 10 18 cm, average 
temperature of the cloud T = 10K, and angular velocity of the cloud £2 = 
10~ 16 rad s _1 . In this case the speed of sound is a s ~ 19200cm s _1 and 
corresponding initial density [8] is p = 10~ 22 g cm~ 3 and accretion rate is 
M = 1.06 x 10 20 g s _1 . In the isothermal phase of the cloud collapse, density 
p oc r~ 2 [9] and the velocity is constant. When opacity becomes high enough 
to trap radiations (say, at r = r tr ), the cloud collapses adiabatically with 
p oc r~ 3 / 2 . In presence of rotation, centrifugal barrier forms at r = r c , 
where centrifugal force balances gravity. Density falls off as p oc r" 1 / 2 in 
this region [10]. The initial constant velocity of infall becomes 8900cm s _1 
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and below r = r c velocity oc r" 1 / 2 was chosen to preserve the accretion rate 
in a disk like structure of constant height. Since for the parameters chosen 
(generic as they are) r c > r tr , we chose T oc 1/r inside the centrifugal 
barrier (r < r c ) as in an adiabatic flow. We follow the collapse till a radius 
of 10 12 cm is reached. 



3. Models and Results 



We use the following models parameterized by a Ad and /. 

Model A: a Ad = l.e - 16 and / = 1. 

Model B: a Ad = l.e - 14 and / = 1. 

Model C: a Ad = l.e - 12 and / = 1. 

Model D: a Ad = 1-e — 10 and / = 1 (same as in Paper I) 

Model E: a Ad = l.e - 16 and / = 100. 

Model F: a Ad = l.e - 14 and / = 10. 

Model G: a Ad = l.e - 12 and / = 5. 



Fig. 1 shows the evolution of adenine abundance X Ad with time (upper 
axis in seconds) and with logarithmic radial distance (in cm). We note that 
generally speaking, at r = 10 16 cm, already the abundance X Ad has reached 
almost the saturated value. This is because, as the upper axis indicates, 
most of the time is spent in this region during collapse. The final abundance 
in different models are: (A) X Ad = 1.36 x 10" 34 , (B) 1.36 x 10" 26 , (C) 
1.36xl0~ 18 , (D) 6.35xl0~ n , (E) 1.2xl0~ 22 , (F) 1.34X10" 20 , and (G) 1.8x 
10~ 14 respectively. Note that when X Ad is really small, it is proportional 
to dr Ad (for fixed /) as expected from a reaction with four sequence (see, 
e.g, Models A, B and C). But when its abundance it significant, HCN 
otherwise participating in other reactions also contribute significantly to 
adenine formation. Similarly, at low X Ad , for / ^ 1, the final abundance is 
proportional to / 6 for the same value of a Ad . 

If the present detectability limit of abundance is around 10~ n [5], then 
it is clear that adenine processed in our method should not be detectable 
(except for Model D) even though it may be enough to contaminate and 
flourish in some planets as we suggest. With a molecular weight of 135 for 
adenine, one could imagine that an abundance of 10 -21 or less should really 
be considered as insignificant as far as the contamination theory goes. In 
that case Models A, B and E must be rejected. This would correspond to 
a lower limit of < a Ad > as ~ 10~ 13 (where we use <> to indicate an 
average over the whole chain of reactions leading the adenine formation 
from HCN.). On the other hand, even when a\ ~< a Ad >, / could be 
large enough to have eventual significant production (Model F). Thus, our 
a — f model implies that one needs to study the reaction rate of not only 
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HCN + HCN -> H2C2N2, but also every stages of HCN addition in order 
to come a definitive conclusion in this regard. 

4. Conclusion 

In presence of radiative association, adenine abundance X^d in an inter- 
stellar cloud seems to be roughly proportional to a\ d f % for small Xj\d- 
This means that the measurements of both a Ad and / must be made very 
accurately. We studied the a — f parameter space and found that while 
some region could produce significant abundance, a smaller region produce 
detectable (with present day technology) amount, while the rest produces 
abundance insignificant enough to dismiss the contamination theory. One 
must wait for the technological advancements to improve laboratory ex- 
periments in extreme conditions and to improve the detectability limit in 
order to come to a firm conclusion. 
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Figure 1.: Evolution of adenine abundance with radial distance and 
time (upper axis) in a a — f model. Various models are marked on the 
curve. See text for details. 



